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ABSTRACT 

We employ Mathematica to find ZAr-invariant subgroups of Eg for application 
in M-theory. These ZTv-invariant subgroups are phenomenologically important 
and in some cases they resemble the gauge groups of our real world. We present 
. a specific example of Z7-invariant subgroups of £"8, which turn up in orbifold 

}^ . compactification of M-theory. However, the procedure can be applied for any 

^ ' group that acts by shifts in the root lattice of semi-simple Lie groups. 

O ■ PACS: 02.20.Rt, 11.25.Mj Keywords: orbifold compactification, 

[ string-theory, M-theory 



p . 1 Introduction 

^ , In models where part of spacetime is compactified, the geometry of compact space affects the gauge 
symmetries of the model. Herein, we consider the Hofava-Witten M-theory [T], where the ll*'^ dimension 
is compactified to an interval, /, and there are two ten-dimensional planes at the boundaries of /. It is 
■ convenient to identify I = /'Ij2, acting as Z2 : — > —(j), so that the boundary of / consists of the fixed 
. points of this Z2-action. On each one of these ten-dimensional spacetime planes, (M^'^)^^^, there is an 
independent copy (a principal vector bundle) of Eg gauge fields. To produce considerably more realistic 
, models with 4-dimensional spacetime, one may proceed as follows: 

o . 

. 1. Impose twisted periodicity conditions on six of the ten dimensions of the boundary spacetime planes, 
passing from C (M"'^'^)^.^ to (T^/A) = ((R^/A)/A), where A is a suitable 6-dimensional lattice 
and A is a symmetry of A. We consider A = Z^v. 
^ I 2. Simultaneously embed the A action into the Eg structure group of the gauge fields on each of the 
I two boundary-spacetimes, the structure groups are broken to subgroups of £"8 that are invariant with 

respect to the A-action. 



This is referred to as "compactifying the Hofava-Witten M-theory on a T^/A orbifold", and A is the 
"orbifold group." Typically, A acts by rotations on the compact space coordinates, and at the same time 
by shifts in the £"8 i^oot lattice. 

In Ref. [2], we have constructed Zy-orbifold models in M-theory. We used Mathematica to find the 
Zy-invariant subgroups of Eg. In this paper we present the details of the Mathematica computation codes 
and the procedure that we have used. This procedure may be used for higher order (iterated) orbifolds as 
well, and also in any situation where one needs to find the Z 7v-invariant subgroups of any of the semi-simple 
Lie group, where Z^r acts by shifts in the root lattice. Since each Lie group factorizes as a product of a 
semi-simple part and its abelian factors, the latter are then treated separately [2]. 
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2 The Algorithm 

Let V denote a Zjv-generating shift vector acting on the root lattice 7' of a group G, where G is one of 
the simple lie groups An, Bn, Gn, Dn, Eq, Ej or E^- The root vectors of G that are invariant with respect 
to the u-projection 

e2™-P |p) = |p) , per, (2.1) 

are the root vectors of a subgroup H C G that is invariant with respect to the whole group generated 
by the element corresponding to v. Different ZTv-generating shift vectors define different subgroups; upon 
identifying those that are equivalent by G-conjugation, we find the inequivalent ZAr-invariant subgroups, 
Hj, with/= 1,2, 

Note: The condition (I2.ip is trivially satisfied for the root vectors = {0,0,- •• ,0} corresponding to 
Cartan generators of G. Therefore, rank(i77) = rank(G), and all Zjv-invariant groups are regular. 

Step 1: Find the set of positive root vectors of G, denoted W. 

Step 2: Based on the restrictions implied, we construct all possible Z^r shift vectors v. 
Step 3: Find all the subgroup^] Hj C G. 

Step 4: For each one of the subgroups Hj C G, define the following four variables: 

t is the set of positive root vectors of Hj C G; 

p := |t| is total number of positive root vectors in Hj C G; 

r := rank'(ff7), defined as the rank of semi-simple part of Hj C G, i.e., without [7(l)-factors; 
m is number of Ai factors, if any, in Hj C G. 

These three variables can be read off by looking at the subgroup Hj and can be used as identifiers of 
the group. If these three variables do not suffice to identify Hj C G unambiguously, define another 
variable: 

7712 is the number of A2 factors in Hj, if any. 

If {p, r, m, turns out not to suffice to identify Hj C G unambiguously, we look for ^43, A^. . . 
factors in Hj, the numbers of which, 771,3, "^4- • • j will be necessary to identify Hj C G unambiguously. 

Step 5: Pick the first v from Step 2. 

Step 5. a: Set t = 0. For all Wa € W, 110 v-Wa = Z, append Wa into the set t. 

Step 5.b: Compute {p, r,m, . . .} of this t (see Section [5] for the procedure). 

Step 5.c: Identify the subgroup Hj C G by comparing {p, r,m, . . .} with the list from Step 4. 

Step 6: Pick the next v from Step 2, and go to Step 5. a. 

Steps 1-4 are preparatory. In particular. Step 4 sets up the string of identifiers {p, r, 77T,, 771,2, ••• } as 
unambiguous addresses of all the subgroups Hj of a given simple Lie group G. In our case, with G = Eg 
and Ziv acting by shifts in the root lattice, {p,r,m} — in many cases only a subset of this — sufficed. 

^For all simple Lie groups of rank < 8 and several of higher rank, the maximal subgroups are listed in the 
literature, such as Ref. [3]. Alternatively, these may be obtained via the technique using extended Dynkin diagrams. 
The non-maximal subgroups are then found by iteratively finding the maximal subgroups of the maximal subgroups, 
and then dropping some of the factors. 

^Since the root lattice shift v corresponds to a generator g{v) e "Zn so that all elements of Zjv are powers oi g{v), 
it follows that root vectors satisfying v wa = Z are in fact invariant with respect to all of Zat. 
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3 Roots and Shift Vectors 



We take the adjoint representation of the group G and calculate it's positive root vectors from the highest 
root Take for example the group G = Eg. The highest root of the irreducible representation 248 of Es 
is {0,0,0,0,0,0,1,0}, which can be read off from the Dynkin diagram of Eg. The entire root system can 
be obtained from the highest root by subtracting from it the positive simple root vectors as follows: in any 
given root vector w, the value of the positive n^^ component indicates how many times the n}^ positive 
root may be subtracted [3j. For example, wi = {2,1,0,0,0,0,0,0} is the first positive root and may be 
subtracted from itself twice, producing (to save space, negative root vector components are denoted by an 

, -r dcf , dof „ , N 

over-bar: 1 = —1, 2 = —2, etc.): 

{2,1,0,0,0,0,0,0} > {0,0,0,0,0,0,0,0} > {2,1,0,0,0,0,0,0}. (3.1) 

All three are indeed in the root system of Eg. Starting with {0,0,0,0,0,0,0,1}, in this way we end with 
{0,0,0,0,0,0,0,1}, having found all 248 roots, spanning 248, the adjoint representation of Eg. 

The procedure actually applies generally to all representations of all simple Lie groups; however, general 
representations are spanned by weight vectors, while the adjoint representation is spanned by roots. By 
plotting the weights (roots) below those from which they are obtained by subtracting positive roots and 
connecting them by arrows, we obtain a "spindle shaped" graph called the weight (root) diagram of the 
representation. In the root diagram, the middle row of the root diagram is populated by eight copies of 
{0,0,0,0,0,0,0,0}, representing the eight Cartan generators. The row immediately above the middle is 
populated by the 8 positive simple root vectors, and the 120 = (248 — 8)/2 roots above the middle row are 
the positive root vectors of Eg. 

The following Mathematica code computes the positive root vectors of Eg. 

Input (1) 

a = {{2, 1,0, 0,0, 0,0,0}, {1,2, 1,0, 0,0, 0,0}, {0,1, 2, 1,0, 0,0,1}, {0,0, 1,2, 1,0, 0,0}, 
{0,0, 0,1, 2, 1,0,0}, {0,0, 0,0, 1,2, 1,0}, {0,0, 0,0, 0,1, 2,0}, {0,0, 1,0, 0,0, 0,2}}; 
g[0] = {{0,0,0,0,0,0,1,0}}; 

g[l] = Table [Flatten [g[0]] - a [ [Flatten [Position [Flatten [g [0] ], 1] ][ [p] ]]] , 

{p , Length [Flatten [Position [Flatten [g [0] ] , 1] ] ] }] ; 
e[x_] := e[x] = Union [Flatten [Table [Table [If [g [x] [ [j] ][ [i] ] == l,g[x][[j]] - a[[i]]], 

{i, 8}],{j, Length [g [x] ]}] , 1]]; 
g[x_] := g[x] = Delete [e[x -1], 1]; 
Flatten [Table [g[m] , {m, 0, 28}], 1] 

Output (1) 

{{0,0, 0,0, 0,0, 1,0}, {0,0, 0,0, 0,1, 1,0}, {0,0, 0,0, 1,1, 0,0}, {0,0, 0,1, 1,0, 0,0}, 
{0,0, 1,1, 0,0, 0,0}, {0,1, 1,0, 0,0, 0,1}, {0,1, 0,0, 0,0, 0,1}, {1,1, 0,0, 0,0, 0,1}, 
{1,0, 0,0, 0,0, 0,1}, {1,1, 1,0, 0,0, 0,1}, {1,0, 1,0, 0,0, 0,1}, {1,0, 1,1, 0,0, 0,0}, 
{1,1, 1,1, 0,0, 0,0}, {1,0, 0,1, 1,0, 0,0}, {1,1, 0,1, 1,0, 0,0}, {0,1, 0,1, 0,0, 0,0}, 
{1,0, 0,0, 1,1, 0,0}, {1,1, 0,0, 1,1, 0,0}, {0,1, 1,1, 1,0, 0,0}, {1,0, 0,0, 0,1, 1,0}, 
{1,1, 0,0, 0,1, 1,0}, {0,1, 1,0, 1,1, 0,0}, {0,0, 1,0, 1,0, 0,1}, {1,0, 0,0, 0,0, 1,0}, 
{1,1, 0,0, 0,0, 1,0}, {0,1, 1,0, 0,1, 1,0}, {0,0, 1,1, 1,1, 0,1}, {0,0, 0,0, 1,0, 0,1}, 
{0,1, 1,0, 0,0, 1,0}, {0,0, 1,1, 0,1, 1,1}, {0,0, 0,1, 0,1, 0,1}, {0,0, 0,1, 1,1, 0,1}, 
{0,0, 1,1, 0,0, 1,1}, {0,0, 0,1, 1,1, 1,1}, {0,0, 0,1, 0,1, 1,1}, {0,0, 1,1, 0,1, 0,1}, 
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To compute all the root vectors of Eq (not just the positive ones), add the Une 

Do [Print [g [m] ] , {m, 0, 28}] 
at the end of Input (1). 



For E^ and Eq the input codes are similar. For Ej, the highest root of the adjoint representation, 133, 
is {1,0,0,0,0,0,0}. Its (133— 7)/2 = 63 positive root vectors of E-j are found from the following code: 

Input (2) 

a = {{2, 1,0, 0,0, 0,0}, {1,2, 1,0, 0,0,0}, {0,1, 2, 1,0, 0,1}, {0,0, 1,2, 1,0,0}, 

{0,0, 0,1, 2, 1,0}, {0,0, 0,0, 1,2,0}, {0,0, 1,0, 0,0, 2}}; 
g[0] = {{1,0,0,0,0,0,0}}; 

g[l] = Table [Flatten [g[0]] - a [ [Flatten [Position [Flatten [g [0] ], 1] ][ [p] ]]] , 

{p , Length [Flatten [Position [Flatten [g [0] ] , 1] ] ] }] ; 
e[x_] := e[x] = Union [Flatten [Table [Table [If [g [x] [ [j] ][ [i] ] == l,g[x][[j]] - a[[i]]], 

{i, 7}],{j, Length [g [x] ]}] , 1]]; 
g[x_] := g[x] = Delete [e[x -1], 1]; 
Flatten [Table [g[m] , {m, 0, 16}], 1] 

For Eq, the highest root (weight of the adjoint representation), 78 is {0 , , , , , 1}. Its (78— 6)/2 = 36 
positive root vectors are found as follows: 

Input (3) 

a = {{2, 1,0, 0,0,0,}, {1,2, 1,0, 0,0}, {0,1, 2, 1,0,1}, {0,0, 1,2, 1,0}, {0,0, 0,1, 2,0}, 

{0,0,1,0,0,2}; 
g[0] = {{0,0,0,0,0,1}}; 

g[l] = Table [Flatten [g[0]] - a [ [Flatten [Position [Flatten [g [0] ], 1] ][ [p] ]]] , 



4 



{p , Length [Flatten [Position [Flatten [g [0] ] , 1] ] ] }] ; 
e[x_] := e[x] = Union [Flatten [Table [Table [If [g [x] [ [j] ][ [i] ] == l,g[x][[j]] - a[[i]]], 

{i, 6}],{j, Length [g [x] ]}] , 1]]; 
g[x_] := g[x] = Delete [e[x -1], 1]; 
Flatten [Table [g[m] , {m, 0, 10}], 1] 

For An, the dimension of the adjoint representation is n(n+2) and the number of positive root vectors 
is {n{n+2) — n)/2 = n(n— 1)/2. The Mathematica code computing the positive root vectors of An, for 
n = 5 for example, is: 

Input (4) 

n = 5; 

d = {{2, -1, 0}, {0, -1, 2}, {-1, 2, -1}}; 

a = If [n > 1, Flatten [{{{PadLeft [d [ [1] ] , n, 0, n - 3]}}, {Table [PadLef t [d [ [3] ] , 

n, 0, n - i - 2] , {i, n - 2}]}, {{PadRight [d [ [2] ] , n, 0, n - 3]}}}, 2], {2}]; 
g[0] = {RotateLeft [PadRight [{1, l}, n, 0], 1]}; 

g[l] = Table [Flatten [g[0]] - a [[Flatten [Position [Flatten [g [0] ] , l]][[p]]]], 

{p, Length [Flatten [Position [Flatten [g [0] ] , 1]]]}]; 
g[x_] := g[x] = Delete [Union [Flatten [Table [Table [If [g[x - l][[j]][[i]] == 1, 

g[x - 1] [[j]] - a[[i]]] , {i, n}] , {j, Length[g[x - 1]]}] , 1]] , 1] ; 
Flatten [Table [g[m] , {m, 0, n - l}] , 1] 

For Bn, the dimension of the adjoint representation is n(2n+l) and the number of positive root vectors 
is (n(2n + 1) — n)/2 = n^. The Mathematica code computing the positive root vectors of Bn, for n = 5 for 
example, is: 

Input (5) 

n = 5; 

d = {{2, -1, 0}, {-1, 2, -2}, {0, -1, 2}, {-1, 2, -1}}; 

a = Flatten [{{{PadLeft [d [ [1] ] , n, 0, n - 3]}}, {Table [PadLef t [d [ [4] ] , n, 0, n - i - 2] , 
{i,n - 3}] }, {{PadRight [d[ [2] ] ,n, 0, n - 3]}}, {{PadRight [d[ [3] ] ,n, 0, n - 3]}}}, 2]; 
g[0] = {PadRight [{0, 1, 0}, n, 0]}; 

g[l] = {Flatten [g[0]] - a [[Flatten [Position [Flatten [g [0] ] , 1]][[1]]]]}; 

e[x_] := e[x] = Union [Flatten [{Table [Table [If [g [x] [ [j] ][ [i] ] == 1, g[x][[j]] - a[[i]], 
If [g[x] [[j]] [[i]] == 2, g[x][[j]] - a[[i]]]], {i, n}], { j , Length [g [x] ]}] , 
Table [Table [If [g[x - l][[l]][[k]] == 2, g[x - 1][[1]] - 2 a[[k]]], {k, n}] , 
{1, Length[g[x - 1]]}]}, 2]] ; 

g[x_] := g[x] = Delete [e[x -1], 1]; 

Flatten [Table [g[m] , {m, 0, 2 n - 2}] , 1] 

Similarly to Bn, the dimension of the adjoint representation of C„ is also n(2n+l) and the number of 
positive root vectors is also (n(2n+l) — n)/2 = v?. The Mathematica code computing the positive root 
vectors of Cn, for n = 5 for example, is: 

Input (6) 

n = 5; 

d = {{2, -1, 0}, {-1, 2, -2}, {0, -1, 2}, {-1, 2, -1}}; 

a = Transpose [Flatten [{{{PadLeft [d [ [1] ] , n, 0, n - 3]}}, {Table [PadLef t [d [ [4] ] , n, 0, 
n - i - 2], {i, n - 3}]}, {{PadRight [d [ [2] ] , n, 0, n - 3]}}, 
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{{PadRight[d[[3]] , n, 0, n - 3]}}}, 2]]; 
g[0] = {PadRight [{2} , n, 0]}; 

g[l] = {Flatten [g[0]] - a [ [Flatten [Position [Flatten [g [0] ] , 2]][[1]]]]}; 

e[x_] := e[x] = Union [Flatten [{Table [Table [If [g [x] [ [j] ][ [i] ] == 1, g[x][[j]] - a[[i]], 
If [g[x] [[j]] [[i]] == 2, g[x][[j]] - a[[i]]]], {i, n}], {j , Length [g [x] ]}] , 
Table [Table [If [g[x - l][[l]][[k]] == 2, g[x - 1][[1]] - 2 a[[k]]], {k, n}] , 
{1, Length[g[x - 1]]}]}, 2]] ; 

g[x_] := g[x] = Delete [e[x -1], 1]; 

Flatten [Table [g[m] , {m, 0, 2 n - 2}] , 1] 

For Dn, the dimension of the adjoint representation is n(2n— 1) and the number of positive root vectors 
is (n(2n— 1) — n)/2 = n(n— 1). The Mathematica code computing the positive root vectors of B„, for n = 5 
for example, is: 

Input (7) 

n = 5; 

d = {{2, -1, 0, 0}, {-1, 2, -1, -1}, {0, -1, 2, 0}, {0, -1, 0, 2}, {-1, 2, -1, O}}; 

a = Flatten [{{{PadLeft [d [ [1] ] , n, 0, n - 4]}}, {Table [PadLef t [d [ [5] ] , n, 0, n - i - 3] , 

{i, n - 4}]}, {{PadRight [d[ [2] ] , n, 0, n - 4]}}, {{PadRight [d [ [3] ] , n, 0, n - 4]}}, 

{{PadRight [d[ [4] ] , n, 0, n - 4]}}}, 2] ; 
g[0] = {PadRight [{0, 1, 0}, n, 0]}; 

g[l] = Table [Flatten [g[0]] - a [[Flatten [Position [Flatten [g [0] ] , l]][[p]]]], 

{p , Length [Flatten [Position [Flatten [g [0] ] , 1] ] ] }] ; 
g[x_] := g[x] = Delete [Union [Flatten [Table [Table [If [g[x - l][[j]][[i]] == 1, 

g[x - 1] [[j]] - a[[i]]] , {i, n}] , {j, Length[g[x - 1]]}] , 1]] , 1] ; 
Flatten [Table [g[m] , {m, 0, 2 n - 4}] , 1] 

The highest root, the level of simple root vectors (i.e., the height of the tower of positive roots) and 
the dimension of the adjoint representation can be found in Table 8 of [3]. Table 9 of Ref. [3] gives the 
positive root systems of a few low-rank simple Lie groups. We leave it to the diligent Reader to adapt the 
above Mathematica codes for the remaining simple Lie groups, G2 and F4. 

— * — 

In constructing T^/Zn = {M.^ / A) /Zj\f orbifolds for superstring theory and its M-theory extension, 
the choices of the Zn shift vectors (representing the embedding in the gauge group) are restricted. For 
example, in M-theory, the shift vectors must satisfy a supersymmetry condition, and in string theory they 
satisfy an additional modular invariance condition. 

We give an example of Z7 vectors. There are 428 eight-component vectors that may be constructed 
with the components taking values in the set {0, ^, |, |, |, |, |}. The supersymmetry restriction requires 
that the components of a Zjv vector add up to an integer. The following is a code that produces all the 
"supersymmetric" Z7 vectors. We have shown only a sample of the output. Note that in order to find all 
the possible vectors preserving supersymmetry, we need to take permutations of each one of the vectors 
shown in this output. 

Input (8) 

a={ 1/7 , 2/7 , 3/7 , 4/7 , 5/7 , 6/7} ; 

b=Flatten [Table [Tuples [a, i] ,{i,2,8}] , 1] ; 

c=Union [Table [Sort [b[[i]]] , {i .Length [b] }] ] ; 
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d=Table [Total [c [ [i] ] ] , {i , Length [c] }] ; 
e=Table [IntegerQ [d[[i]]] , {i .Length [d] }] ; 
f=Flatten [Position [e , True] ] ; 
g=Table [c [ [f [ [i] ] ] ] , {i , Length [f ] }] ; 
q=Table [PadRight [g [ [i] ] , 8] , {i , Length [g] }] 
"Total no. of Z7 Vectors" 
Length [q] 

Output (8) 

{{i,f ,0,0,0,0,0,0}, {f, 1,0,0,0,0,0,0}, {f, 1,0,0,0,0,0,0}, {i,i,f ,0,0,0,0,0}, 
{i,f , 1,0,0,0,0, 0}, {i,f,f ,0,0,0,0,0}, {f,f,f ,0,0,0,0,0}, {f,f,f ,0,0,0,0,0}, 
{f,f,f ,0,0,0,0,0}, {|,|,f ,0,0,0,0,0}, {|,f,f ,0,0,0,0,0}, 1,0,0,0,0}, 



r3 366666 61 r3 444446 6l f3 444455 61 r3 444555 51 

^T'T'T'T'T'T'T'T-'' L7''7''7''7'7'7''7'7'J' \7'7'7''7''7''7'7''7'J' L'7''7'7''7'7'7''7''7/' 

/3 456666 6i r3 555666 6i r4 444445 6i ri 444455 5; 

'^7'7'7'7'7'7'7'7J' XysyjyjyjyJY'S'j'jy/j Xysyjyjyjyj'j'jyjyjj Xysyjyjyjyjyj'j'j'j'/) 

r4 446666 6i /i 455666 61 /i 555566 6\ f5 555556 6i\ 

l.7'7'7'7'7'7'7'7J' 17'7>7'7>7'7'7>7J> 17'7>7'7>7'7>7'7J> 17'7>7'7>7'7>7'7JJ 

"Total No. of Z7 Vectors" 
428 



One may use a similar code for generating general Z^r shift vectors in the root lattice. 



4 Subgroups of G 

Our next step is to find all the regular subgroups of G. This can be done manually, for example, by 
using the extended Dynkin diagram technique [1]. The procedure is to remove nodes from the extended 
Dynkin diagram in all possible ways, and read off the subgroup represented by the remainder. For example 



Qr^V^V^V^V-Or-c^7 o„ o— o— 6— 

a'^ a a° a" a' a" a"^ a a° a° 

Figure 1: Nodes , o? and removed from the extended Dynkin diagram of gives the regular 
subgroup 5*010 S\J2 x \J\\ the number of C/(l) factors is determined from the difference in the rank. 



(see Fig [T|), if we take out the nodes a^, c? and from the extended Dynkin diagram of E'g, we get 
SO\Q X SU2 X as a subgroup. Following this procedure, we get a total of 99 subgroups of £^8 given in 
Tabled! 

Alternatively, we can use Mathematica to find all these subgroups of Eg: 
Input (9) 

b = Flatten [{Table [a [i] , {i, 8}], Table [d[j], {j, 4, 8}] ,Table [e [k] , {k, 8}]}]; 
c = Flatten [Table [Tuples [b, i] , {i, 8}], 1]; 
f = Union[Table[Sort[c[[i]]] , {i. Length [c] }]] ; 

g = Table [Sum [f [[j]] [[i]] [[1]] , {i, 1, Length [f [ [j] ]]}], {j , Length[f]}]; 
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h = Table [If [g[[i]] <= 8, True, False], {i. Length [g] }] ; 

m = Flatten [Position [h, True]]; 

s = Table[f [[m[[i]]]] , {i, Length [m] }] 

In the above code, we look for possible combinations of An < n < 8), Dn (4 < n < 8) and En 
(6 < n < 8). This code may not work for computers low in processor power and memory. In such a 
case, one might want to break down the evaluation in k pieces (k is equal to the rank of the group G) by 
changing the variable i in the second line of the input as {i,k,k}, instead of {i,8} and collecting all the 
results at the end. 

Using the Dynkin Diagram technique or Mathematica, one can find all the regular subgroups of G in 
general. 

Table 1: All the subgroups of Eg and the values of the identifiers. 





Subgroups 


p 


m 


r 


c 




Subgroups 


p 


m 


r 


c 


1 


Es 


120 





8 




51 


SOs X sui 


14 


2 


6 


/ 


2 


Ej X SU2 


64 


1 


8 




52 


SU5 X SUs X SU2 


14 


1 


7 


/ 


3 


Er 


63 





7 


/ 


53 


SU5 X SC/| 


14 


4 


8 


/ 


4 


SOlQ 


56 





8 




54 


sui X sui 


14 


2 


8 


/ 


5 


SOu X SU2 


43 


1 


8 




55 


SOs X SU2 


13 


1 


5 




6 


50i4 


42 





7 


/ 


56 


SU5 X SU3 


13 





6 




7 


Eg X SU3 


39 





8 




57 


SU5 X SU'I 


13 


3 


7 




8 


Eg X SU^ 


38 


2 


8 




58 


sui X SU2 


13 


1 


7 




9 


Eq X SU2 


37 


1 


7 


/ 


59 


SUi X SUl X SU2 


13 


1 


8 




10 


E(> 


36 





6 


/ 


60 


SOs 


12 





4 




11 


SUg 


36 





8 


/ 


61 


SU5 X sui 


12 


2 


6 




12 


SO12 X SUs 


33 





8 




62 


sui 


12 





6 




13 


SO12 X sui 


32 


2 


8 




63 


sUi X sui 


12 





7 




14 


SO12 X SU2 


31 


1 


7 




64 


sUi X SU3 X sui 


12 


3 


8 




15 


SO12 


30 





6 


/ 


65 


sui 


12 





8 




16 


SUs X SU2 


29 


1 


8 




66 


SU5 X SU2 


11 


1 


5 




17 


SUs 


28 





7 


/ 


67 


SU4 X SU3 X sui 


11 


2 


7 




18 


SOio X SU4 


26 





8 




68 


SU4 X sui 


11 


5 


8 




19 


SOio X SU3 X SU2 


24 


1 


8 




69 


sui X sui 


11 


2 


8 




20 


SUj X SUs 


24 





8 




70 


SU5 


10 





4 




21 


SOw X SUi 


23 





7 


/ 


71 


SUi X 5C/3 X SU2 


10 


1 


6 




22 


SOio X SUl 


23 


3 


8 


/ 


72 


SUi X sui 


10 


4 


7 




23 


SU7 X 5C/| 


23 


2 


8 


/ 


73 


sui X SU2 


10 


1 


7 




24 


SOio X sui 


22 


2 


7 


/ 


74 


sui X sui 


10 


4 


8 




25 


SOs X 5;75 


22 





8 


/ 


75 


SU4 X SU3 


9 





5 





Continued. . . 
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Table 1: All the subgroups of Eg and the values of the identifiers. (Continued) 





Sn heTonns 




ni 


I" 


c 




Sn heTonns 




ni 


r 




26 


sr/v X suo 


22 


1 


7 


/ 


76 


SUa X SU^ 


9 


3 


6 




27 




21 


1 


6 


/ 


77 




9 





6 




28 




21 





6 


/ 


78 


SU? X 


9 


3 


7 




29 


.ST/r X SUa 


21 





8 


/ 


79 


sU'i X 5r/a 


9 


6 


8 




30 


.Win 


20 





5 




80 


5r//i X 5r/o 


8 


2 


5 




31 




20 





8 




81 


5772 X SUn 


8 


2 


6 




32 


SOn X Sf/zi X Sf/o 


19 


1 


8 




82 


SC/t X sul 


8 


5 


7 




33 


SL^fi X X SUo 


19 


1 


8 




83 




8 


8 


8 




34 


SOn X S'[//i 

X x\ i_/ ij. 


18 





7 




84 


SUa X Sf/o 


7 


1 


4 




35 


.Wa X .ST/n^ 


18 





8 




85 


sul X SC/o 


7 


1 


5 




36 


SC/fi X SC/ci 


18 





7 




86 


kJ \-J ^ W 


7 


4 


6 




37 


Sf/fi X sr/i? 


18 


3 


8 




87 




7 


7 


7 




38 


50s X X SU'^ 


17 


2 


8 




88 


5'C/4 


6 





3 




39 


5C/fi X sul 


17 


2 


7 




89 


sul; 


6 





4 




40 


iSTA X SUa X 5r/o 


17 


1 


8 




90 


.Qr/o X .ST/a 


6 


3 


5 




41 


.Wa X 5r/q X .ST/o 


16 


1 


7 


/ 


91 


1JU2 


6 


6 


6 




42 


SOs X SUn 


16 


4 


8 


/ 


92 


SU-x X Sf/o 


5 


2 


4 




43 


SUk X ST/o 


16 


1 


6 


/ 


93 




5 


5 


5 




44 


SU^ X SfJ/i 


16 





7 


/ 


94 


SfA X SUo 


4 


1 


3 




45 


SU^ X SC/o^ 


16 





8 


/ 


95 


SU^ 


4 


4 


4 




46 


SOfi X S'C/q 


15 





6 


/ 


96 


SU-x 

KJ ^ 


3 





2 




47 


SOs X 


15 


3 


7 


/ 


97 


sul 


3 


3 


3 




48 




15 





5 


/ 


98 


sul 


2 


2 


2 




49 


X SU3 X 5C/| 


15 


2 


8 


/ 


99 


SU2 


1 


1 


1 




50 


5[7| X SU3 


15 





8 


/ 


100 
















The C/(l) factors are not shown; regular subgroups have (8— r) of them. 



Once we get all the subgroups of E'g, we calculated the number of positive root vectors for each subgroup 
and list them in column p. The values of the other identifiers (r, m and possibly m2, ms, . . . ) turned out 
not to be necessary in most cases for our purpose^: Before using them, we found the possible candidates 
which are Z7 invariant subgroups of £"8 through a procedure given in Input/Output (10). This greatly 
reduced the complexity of the codes in the next section and saves in the Mathematica evaluation time. 

We have 428 Z7 shift vectors in Output |^ and once we take their permutations, this gives a total 
of 823,542 shift vectors. We take the first Z7 vector {^,1, 0,0, 0,0, 0,0} from the previous section and 

■'The identifiers r and m are shown in Table [T] for completeness, and for the benefit of possible generalizations to 
Zjv-actions where the supersymmetry condition is relaxed. The additional identifiers, mi in Step 4, are easily added. 
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calculate the number of positive root vectors that satisfy the condition p v = Z using the following code: 
Input (10) 

q = [not shown here: 4^8 Z7 vectors from Output ^B^ ) 

w = [not shown here: 120 positive root vectors of from Output (Cp ) 

p = FlattenCTable [Permutations [q[[i]]] , {i, 1, l}] , 1]; 

a = Length [w] ; 

1 = Length [p] ; 

u = Table [Table [w [ [i] ] .p[[j]] , {i, a}], { j , l}] ; 

V = Table [Table [IntegerQ [u [ [j , i]]], {i, a}], { j , l}] ; 

r = Table [Count [v [ [j] ] , True], { j , l}] ; 

Union [r] >>>Z7_Roots; 

Output (10) 
{37,42} 

The Output of this evaluation ({37,42}) is written in an external file ZVJioots. We do this evaluation 
for the other Z7 vectors in Output (B^ and the results are collected from the text file ZVJioots. This gives 
the possible values of p for Z7 vectors as 



This narrows down our choices to 30 subgroups of Eg, marked by a check in column c of Table [H Now we 
use the values of m (number of Ai factors in Hj C Eg) to identify the possible subgroups Hj C Eg. When 
p and m do not specify Hi C Eg unambiguously, we use the values of r (rank of the semi-simple part of 
Hi). The values of p, m and r are also calculated from the root vectors that survive a Z7 shift. This is 
shown in the next section. 

5 Zjv Invariant Subgroups of G 

To illustrate the procedure of calculating the values of m and r from the surviving root vectors we give 
the same example as in Section 3 of our previous paper [2]. Take the shift vector v = {i, i, 0, |, 0, 0, |, 0}, 
which is one of the permutations of {4, = ,0,0,0,0}. The E^ root vectors that survive the shift are 



and are thus invariant under the action of the whole Z7 generated by this shift. Call these root vectors 
i = 1, 2, • • • 15, and set p = 15. 

Next, we need to identify which subgroup Hi C E^ — from among those listed in Table [1] — do these 
roots (together with their negatives and the Cartan root vectors) generate. We look for possible relations 



p : {14, 15, 16, 21, 22, 23, 28, 30, 36, 37, 42, 63} 



(4.1) 



{0,0,0,0,1,1,0,0}, {1,1,0,0,0,0,0,1}, {1,1,1,0,0,0,0,1}, 

{1,1,0,0,1,1,0,0}, {0,0,1,0,1,0,0,1}, {0,0,0,0,1,0,0,1}, 

{1,1,0,0,0,1,0,0}, {0,0,0,0,0,1,0,1}, {0,0,1,0,0,1,0,1}, 

{0,1,1,1,1,1,1,0}, {1,0,1,1,0,0,1,0}, {1,1,1,0,1,0,0,0}, 

{1,1,1,0,1,1,0,1}, {1,1,0,0,1,1,0,1}, {0,0,1,0,0,0,0,2}, 



(5.1) 
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in the form t[i] + t[j] = t[k] and find the following: 



t[2] + t[U] = t[l] , t[3] + t[l3] = t[l] , t[3] + t[15] = t[2] , (5.2a) 

t[5] + t[9] = t[l] , t[5] + t[l2] = t[2] , t[6] + t[8] = t[l] , (5.2b) 

t[6] + t[12] = t[3] , t[6] + t[l5] = t[5] , t[7] + t[8] = t[2] , (5.2c) 

t[7] + t[9] = t[3] , t[7] + t[13] = t[5] , t[7] + t[14] = t[6] , (5.2d) 

t[9] + 1[15] = t[8] , t[10] + t[n] = t[4] , t[12] + 1[13] = t[8] , (5.2e) 

t[12] + t[U] = t[9] , t[U] + t[15] = t[13]. (5.2f) 



Note that the root vectors t[7], t[10], t[12], t[14] and t[15] can not be expressed as a sum of any 
other root vectors, which means that they must correspond to 6 positive, simple root vectors in Hj, whence 
the rank of the semi-simple part of Hj must be 6, and the remaining two Cartan roots correspond to a 
U{1)'^ factor. Also, all the 15 root vectors appear in (|5.2|) . meaning there are no Ai factors in Hj, each 
of which would have had to have a single, isolated, positive root vector. So looking at the surviving root 
vectors we can define the variables m and r as: 

m is the number of root vectors that do not appear in the equation of the form t[i] + t[j] = t[k] 
and so must be single, isolated, positive root vectors; here, m = 0. 

r is the number of root vectors that do not appear on the right side of the equation of the form 
t[i\ + t[j] = t[k] and so must be simple; here, r = 6. 

Using {p, m,r} = {15,0,6}, we identify unambiguously the subgroup from Table [1] as SO^ x SU^. 

We employ this analysis in the construction of the Mathematica codes below and find the subgroups 
of £"8 that are invariant under a Zy shift listed in Table [21 and so the whole Z7 group action generated by 
that shift. This reduces our choice to only 14 subgroups. 



Table 2: Zy-invariant subroups of Eg- 





Group 




Group 




Group 




Group 


1 


Er 


5 


SOu 


9 


SUs 


13 


SU5 X SU4 


2 


Ee X SU2 


6 


SOw X SU3 


10 


SUj X SU2 


14 


SU5 X SU3 X SU2 


3 


Eg 


7 


SOio X SU2 


11 


SU7 






4 


SOu 


8 


SOs X SU3 


12 


SUe X SU2 







Input (11) 

q = (not shown here: J^28 Z7 veetors from Output (0) ) 

w = {not shown here: 120 positive root vectors of Eg from Output ^^ ) 

CleanSlate [] ; 

p = Flatten [Table [Permutations [q [ [i] ]] ,{i, 1,1}] ,1] ; 

u = Table [Table [w[ [i] ] .p[[j]] ,{i,Length[w]}] ,{j ,Length[p]}] ; 

V = Table [Table [IntegerQ [u [ [ j , i] ] ] , {i , Length [w] }] , { j , Length [p] }] ; 

s = Table [Flatten [Position [v [ [k] ], True] ], {k, Length [v] }] ; 

t = Table [Table [w[[s[[j]] [ [i] ] ] ] , {i ,Length[s [ [j] ] ] }] , {j ,Length[v] }] ; 
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$ [k_] : =$ [k] =Evaluate [b=Table [Table [t [ [k] ] [ [i] ] +t [ [k] ] [ [ j ] ] , 

{i , Length [t [ [k] ] ] }] , { j , Length [t [ [k] ] ] }] ; 
c = Table [Table [MemberQ [t [ [k] ] ,b[[i,j]]] , {i .Length [t [ [k] ]] }] , {j , Length [t [ [k] ]] }] ; 
f = Position [c , True] ; 

g = Union[Table [Sort [f [[i]]] ,{i, Length [f]}]] ; 
X = Table [g [ [i] ] [[1]] ,{i, Length [g]}] ; 
y = Table [g [ [i] ] [[2]] ,{i, Length [g]}] ; 

h = Table[t[[k]] [[x[[i]]]]+t[[k]] [ [y [ [i] ]]],{i, Length [x]}] ; 

z = Flatten[Table[Position[t[[k]] ,h[[i]]] ,{i,Length[h]}]] ; 

o = Table [1,{1, Length [t [ [k] ]]}] ; 

m = Length [Complement [o, Union [x,y,z]]] ; 

n = Length [Complement [o , z] ] ; ] ; 

Table [If [Length [t [ [k] ] ] ==14 , Evaluate [$ [k] ; If [m==l , a [1] a [2] a [4] , If [m==2 , 
If [n==6 , a [1] a [1] d [4] , If [n==8 , a [1] a [1] a [3] a [3] ] ] , If [m==4 , a [1] a [1] a [1] a [1] a [4] ] ] ] ] , 
If [Length [t [ [k] ] ] ==15 , Evaluate [<I> [k] ; If [m==0 , If [n==5 , a [5] , If [n==6 , a [2] d [4] , 
If [n==8,a[2]a[3]a[3]]]] ,If [m==2, a [1] a[l] a [2] a[4] , If [m==3,a[l] a [1] a[l] d[4] ] ] ] ] , 
If [Length [t [ [k] ] ] ==16 , Evaluate [$ [k] ; If [m==0 , If [n==7 , a [3] a [4] , If [n==8 , a [2] a [2] a [4] ] ] , 
If [m==l , If [n==6 , a [1] a [5] , If [n==7 , a [1] a [2] d [4] ] ] , If [m==4 , a [1] a [1] a [1] a [1] d [4] ] ] ] ] , 
If [Length [t [ [k] ] ] ==21 , Evaluate [$ [k] ; 
If [m==0 , If [n==6 , a [6] , If [n==8 , a [3] a [5] ] ] , If [m==l , a [1] d [5] ] ] ] , 
If [Length [t [ [k] ] ] ==22 , Evaluate [$ [k] ; 
If [m==0 , a [4] d [4] , If [m==l , a [1] a [6] , If [m==2 , a [1] a [1] d [5] ] ] ] ] , 
If [Length [t [ [k] ] ] ==23 , Evaluate [^> [k] ; 
If [m==0 , a [2] d [5] , If [m==2 , a [1] a [1] a [6] , If [m==3 , a [1] a [1] a [1] d [5] ] ] ] ] , 
If [Length [t [ [k] ] ] ==28 , a [7] , If [Length [t [ [k] ] ] ==30 , d [6] , 
If [Length [t [ [k] ] ] ==36 , Evaluate [^> [k] ; 
If [n==6 , e [6] , If [n==8 , a [8] ] ] ] , 
If [Length [t [ [k] ] ] ==37 , a [1] e [6] , If [Length [t [ [k] ] ] ==42 , d [7] , 
If [Length [t [ [k] ] ] ==63 , e [7] ]]]]]]]]]]]], {k , Length [t] }] ; 
Union [%] »>Z7_Groups ; 

Output (11) 
{d[7], a[l]e[6]} 

In the above code Length [t [ [k] ] ] is the number of surviving root vectors (p), m is the number of Ai 
factors (m) and n is the rank (r) of a group. These three variables are calculated from the surviving root 
vectors as explained in the example, £^8 — ^ 50(8) x SU{3) branching. The output of this evaluation is 
written in an external file Z7_Groups where a group An is identified as a[n] , Dn as d[n] and as e [n] . 

For other orbifolds there are situations where p, m and r do not suffice to specify the group unam- 
biguously. In those cases we we look for A2, A3, ■ ■ ■ factors in Hj by looking at root vector relations. For 
example, an A2 factor would have to be spanned by three root vectors that satisfy an equa- 

tion of the form t[i] + t[j] = t[k] and occur in no equation involving any other root vectors. Equivalently, 
we can look for root vectors that do not appear in any equation of the form t[i] + t[j] + t[k] = t[l]. 
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6 Automation 



Due to limitations of computer's processor speed and memory, it may be necessary to partition the com- 
putation. The following shows how it may be done for the TL-j orbifold example in M-theory. 

(i) Collect all the Z7 vectors q in Output ^ and all the Eg root vectors w in Output (OP of section [3] and put 
them in a notebook, say, NB_0. Use the package 'CleanSlatelf] and put this in one of Mathematical s home 
directory. This package helps in clearing the Mathematica kernel memory so that successive evaluations 
can use the maximum possible memory. The input of NB_0 are as follows: 

Input (12) 

q = ; (no output shown here: ^7 vectors from Output ^ ) 

w = ; (no output shown here: 120 positive root vectors of Eg from Output ) 

« CleanSlate .m; 

orbifold = EvaluationNotebook [] ; 
NotebookSave [orbifold] 
NotebookOpen["NB l.nb"] 

(ii) We create a notebook Z7_Generic which contains the code of Input HO]) with some added lines of 
codes to make use of the automation process: 

Input (13) 

NotebookClose [orbifold] 
CleanSlate [] ; 

p = Flatten [Table [Permutations [q[ [i] ]] , {i, a, a}], 1]; 
a = Length [w] ; 
1 = Length [p] ; 

u = Table [Table [w[ [i] ] .p[[j]] , {i, a}], { j , l}] ; 

V = Table [Table [IntegerQ [u [ [j , i]]], {i, a}], { j , l}] ; 

r = Table [Count [v [ [j] ] , True], {j , l}] ; 

Union [r] »> Z7_Roots; 

orbifold = EvaluationNotebook [] ; 

NotebookSave [orbifold] 

7 = a + 1 ; 

"NB " ToString[7] <> ".nb"; 
InputForm [%] 
NotebookOpen [%] ; 

Next we create a notebook Z7_Generator with the following set of codes, 
Input (14) 

Do [NotebookPut [NotebookGet [First [Notebooks ["Z7 Generic .nb"] ] ] / . "«"->/?] ; 
NotebookSave [SelectedNotebook[] ,"NB "oToString [/?]<>" .nb"] ; 
Pause [2] ; 

NotebookClose [SelectedNotebook [] ] , /3 , 1 , 428] 

Once the Input is run, it creates 428 notebooks with the contents of Input ^3\) where the value of 
a = 1, 2, 3, • • • , 428, respectively, for each notebook. The files are created in the $HomeDirectory. 

"^This package is available on-line at: http://library.wolfram.com/infocenter/MathSource/4718/. 
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(iii) Our next step is to evaluate these 428 notebooks in a way such that when we open NB_0^ it au- 
tomatically evaluates it's content and the contents of notebooks NB_1^ NB_2 and so on so forth. The 
NotebookClose [orbif old] input line closes the previous notebook that has been evaluated. In this way 
the screen is not cluttered with open Mathematica notebooks, improving the performance of the computer's 
memory. The memory is also managed by the input line CleanSlate [] . Note that the 'CleanSlate' package 
is called in after Mathematica stores the values of q and w in its memory which is necessary for the whole 
evaluation process. The end result is collected from the text file ZVJioots created in $HomeDirectory and 
is given in Eq. (|4.ip . 

(iv) We apply a similar procedure for the evaluation of the Z7 invariant groups, Input The results 
are collected from the text file Z7_Groups and are summerized in Table [2l 

In order for the automation process to work we need to make the following changes to Mathematica 
preferences, 

1. Notebook Options — ?> File Options — )■ Notebook Autosave (False — ?> True) 

2. Notebook Options — > File Options — )■ ClosingAutosave (False True) 

3. Notebook Options — ?> File Options — )• AutogeneratedPackage (Manual None) 

4. Notebook Options — > Evaluation Options — > Initialization CellEvaluation (Automatic — )• True) 

5. Notebook Options Evaluation Options Initialization CellWarning (True — > False) 

6. Cell Options — >■ Evaluation Options — > Initialization Cell (False True) 

This automation process has been tested in Mathematica version 5.2, 6.0.1 and 7.0. There appears to 
be a problem in the latter two versions, which prevents the evaluation of a notebook when it is opened 
by another notebook, even though the Initialization CellEvaluation and Initialization Cell are changed to 
True (globally). In those versions of Mathematica, the automation process (iii) can be performed using a 
code such as: 

Input (15) 

nb = NotebookOpen ["notebook. nb"] ; 
SelectionMove [nb. All, Notebook]; 
Select ionEvaluate [nb] ; 
orbif old = EvaluationNotebook[] ; 
NotebookSave [orbif old] ; 
NotebookClose [orbif old] ; 

Corresponding changes need to be made also in Input [77]] and Input ( f73|) for the automation process to 
work. 

7 Conclusion 

We have shown in detail how to find the Z7-invariant subgroups of using Mathematica. These groups, 
obtained in orbifold M-theory, turn out to be closely related to string theory compactification down to four 
dimensions: In the limit x^^ 0, the two Zy-invariant subgroups of E^ (one on each of the two boundaries 
of x^^) coalesce into Hj^i^ x Hji ^, which turn out to coincide with the gauge groups found in Zy-orbifold 
models in string theory [S]. We have tested our codes also for Z2, Z3, Z4 and Zg orbifolds. The so-obtained 
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subgroups upon the limit x — )> coincide with those found in string theory compactification. This would 
imply that our codes can also be used for Zg and Z12 orbifolds. 

In the presence of gauge background fields (Wilson lines) the four-dimensional gauge group breaks 
down to some smaller groups. Since these Wilson lines provide additional shifts in the group lattice, it 
should be possible to employ our procedure also in those types of models. 

For the simple Lie groups A„, S„, C„, D„, £^6 and Ej, our procedure can be applied in finding the 
unbroken gauge symmetry under any Zjy rotations. In section [3l we provided the root vectors for these 
groups. As semi-simple Lie groups are products of simple Lie groups, the procedure presented herein 
merely needs to be applied to each factor separately. Finally, the remaining abelian factors in a Lie group 
are then treated separately [2] . 

Finally, our present goal was the demonstration that Mathematica can be used to compute A-invariant 
subgroups of semi-simple Lie groups. Having been motivated by applications in M-theory, we have re- 
stricted ourselves to "supersymmetric" A-actions, and for simplicity to A = Zjv. Generalizations in both 
respects seem to be desirable, but are beyond our present scope. Similarly, it would seem desirable to re- 
structure, optimize and package the computations presented herein into a single, user-friendly interactive 
Mathematica package, but that too is deferred to a subsequent effort. 
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